arXiv: 1503.04518vl [hep-th] 16 Mar 2015 


Noname manuscript No. 

(will be inserted by the editor) 


Effect of phantom dark energy on the holographic 
thermalization 


Xiao-Xiong Zeng^i^, Xin-Yun Li-Fang Li (corresponding author)^ 

^ School of Science, Chongqing Jiaotong University, Chongqing, 400074, China 

^State Key Laboratory of Space Weather, Center for Space Science and Applied Research, Chinese Academy of Sciences, 
Beijing, 100190, China 


the date of receipt and acceptance should be inserted later 

Abstract Gravitational collapse of a shell of charged 
dust surrounded by the phantom dark energy is probed 
by the minimal area surface, which is dual to probe the 
thermalization in the boundary quantum field by ex¬ 
pectation values of Wilson loop in the framework of the 
AdS/CFT correspondence. We investigated mainly the 
effect of the phantom dark energy parameter and chem¬ 
ical potential on the thermalization. The result shows 
that the smaller the phantom dark energy parameter 
is, the easier the plasma thermalizes as the chemical 
potential is fixed, and the larger the chemical potential 
is, the harder the plasma thermalizes as the dark en¬ 
ergy parameter is fixed. We get the fitting function of 
the thermalization curve and with it, the thermalization 
velocity and thermalization acceleration are discussed. 

Keywords Holographic thermalization, phantom 
dark energy, non-equilibrium 

1 Introduction 

Recent years, more and more theoretical physicists have 
paid their attention on the applications of the AdS/CFT 
correspondence [T]. It can not only check the effective¬ 
ness of this correspondence indirectly, but also provide 
a method to deal with some problems in strongly cou¬ 
pled systemjSlEllllElinilZ]- One of the strongly coupled 
systems is the quark gluon plasma which are produced 
in heavy ion colliders such as the RHIC and LHC. The 
properties of the quark gluon plasma have been inves¬ 
tigated extensively and we know that it behaves as an 
ideal fluid with a very small shear viscosity over entropy 
density ratio [S]. But the process of formation of quark 
gluon plasma after a heavy ion collision, often referred 
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to as thermalization, is not well understood until now. 
Recently some physicists have devoted themself to ad¬ 
dressing this problem from the viewpoint of holography, 
which is called holographic thermalization. The most 
prominent property for the thermalization is that it is 
a non-equilibrium process. From the AdS/CFT corre¬ 
spondence, we know that the initial state before the 
thermalization is dual to a pure AdS, and the last state 
after the thermalization is dual to a stationary black 
hole. Therefore, to describe the thermalization process, 
one should construct a dynamical background in the 
bulk, which can be described as black hole formation 
or black hole merger. 

Now, there have been many models to study this far- 
from-equilibrium thermalization behaviors p i 10 lITTlIT^ 
[uniiisiiiain]. A slightly different and simpler model 
on this topic was put forward by Balasubramanian et 
al. [TSlfT^ , where the dynamical background was treated 
as the gravitational collapse of a shell of dust. They 
claimed that the thermalization process can be probed 
by the equal-time two-point correlation functions of 
local gauge invariant operators, expectation values of 
Wilson loop operators, and entanglement entropy, which 
in the gravity side correspond to minimal lengths, ar¬ 
eas, and volumes in AdS space, respectively. An im¬ 
portant conclusion from their model is that the ther¬ 
malization time is closer to the experiment data pro¬ 
duced in RHIC and LHC. In addition, they found that 
the thermalization is a top-down process, in contrast to 
the predictions of bottom-up thermalization from per¬ 
turbative approaches dO], and there is a slight delay in 
the onset of thermalization. Now such model has been 
generalized to the bulk geometry with electrostatic po¬ 
tential high curvature corrections 

[77] , and some other gravity models ^51(^1^13 1 113211551 
[551i551[5in[571l551l551l5?5l5Tll57l551li^ . 
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In this paper, we intend to study holographic ther- 
malization in the bulk surrounded by phantom dark en¬ 
ergy. The phantom dark energy is the scalar field with a 
negative kinetic term with an equation of state param¬ 
eter a; < —1. As already known, the presence of cos¬ 
mological phantom fields continues to receive supports 
from both collected observational data [IS] and theoret¬ 
ical models [IS] . All these programs have pointed out an 
accelerated expansion of the universe, dominated by an 
exotic fluid of negative pressure. Furthermore, there are 
evidences suggesting the exotic fluid could be of phan¬ 
tom nature mM- Since then, an interest in phantom 
fields has grown and resulted in many phantom black 
hole solutions [4^152] . In recent years, many issues per¬ 
taining to phantom black hole, such as thermodynamic 
stability [SS] and light paths [51I5S] have been dealt 
with. In this paper, we study holographic thermaliza- 
tion in the spacetime dominated by phantom dark en¬ 
ergy. The authors in ref. m has studied the holographic 
thermalization dominated by quintessence dark energy 
and their results show that the smaller the state pa¬ 
rameter of quintessence is, the harder the plasma to 
thermalize. Considering the first year WMAP data to¬ 
gether with the 2dF galaxy survey and the supernova 
la data favor the phantom energy over the cosmological 
constant and the quintessence, here we study the effect 
of phantom dark energy on the holographic thermal¬ 
ization. Our results show that the sate parameters of 
dark energy play different roles in the process of holo¬ 
graphic thermalization. Our results here shows that the 
smaller the phantom dark energy parameter is, the eas¬ 
ier the plasma to thermalize when the chemical poten- 
tieal is fixed, which is different from the case produced 
by quintessence dark energy. 

The structure of this paper is outlined as follows. In 
Section we will briefly review the gravitational col¬ 
lapse solution in the spacetime dominated by the phan¬ 
tom dark energy. In Sectior(31 we derive the equations of 
motion of the minimal area surface theoretically firstly, 
and then perform a systematic analysis about how the 
dark energy parameters and chemical potential affect 
the thermalization time resorting to numerical calcu¬ 
lation. We also study the thermalization velocity and 
thermalization acceleration with the fitting functions of 
the thermalization probes. The last section is devoted 
to discussions and conclusions. 


2 Gravitational collapse surrounded by 
phantom dark energy 

Start from the action 

5 = j d^x^[R + 27jF^^'^F^,+ 2A], (1) 


reference recently obtained a phantom Reissner- 
Nordstrom AdS black hole solution 


ds^ =—f{r)dt^ + f ^ (r)dr'^ Fr'^dcf)'^+r^ siii^ (f)dip^, (2) 


where 


2M 


2 2 
a r 
f{r) = 1- + 


( 3 ) 


In which L is the radius of the AdS black hole, rj is the 
phantom dark energy parameter. For ry = 1, the solu¬ 
tion in Eq. represents for the Reissner-Nordstrom- 
AdS black hole. There is also a purely electric gauge 
field given by 


At = {u-^), 
r 


( 4 ) 


where m is a constant corresponding to the electrostatic 
potential at r —>■ oo, which will be related to the chem¬ 
ical potential in the dual gauge theory according to the 
AdS/CFT correspondence. It is defined such that the 
gauge field vanishes at the horizon, namely u = q/vh- 
In the context of the AdS/CFT correspondence, one 
is often interested in dual gauge theories living on flat 
space. In order to achieve this goal, one often use the 
infinite volume limit m- After this step, we can get the 
phantom Reissner-Nordstrom-AdS black brane solution 


ds^ = -F{r)df -t F-\r)dr^ -t r^dx^, (5) 


in which i = 1,2, and 


F{r) = 




( 6 ) 


In this case, the event horizon, defined by F(rh) = 0, is 
plane instead of sphere. According to the definition of 
the surface gravity, we also can get the Hawking tem¬ 
perature emitting from the black brane 


K I 2M 2q2 2rh. 


( 7 ) 


which in the framework of AdS/CFT can be viewed 
as the equilibrium temperature of the dual field theory 
living on the boundary. 

To get a Vaidya type evolving black brane, we will 
reexpress the metric in Eq.© in Eddington-Einkelstein 
coordinate dv = dt + and then rewrite it with 

the inverse radial coordinate z = —. The background 
spacetime in Eq.® changes into 

ds^ = — \—H{z)dv^ — 2dz dv + dxj\ , (8) 

where L has been set to 1, and 


H{z) = 1 - 2Mz^ + 


(9) 
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Treating the mass parameter in Eq.® as an arbitrary 
function of v, Eq.® can be regarded as a gravita¬ 
tional collapse solution surrounded by phantom dark 
energy [I^[T^I21 II 22] , This solution is not a solution of 
the action in Eq. 0 anymore. It inclndes the contribu¬ 
tions of some external matter fields. As one can show, 
such a metric is sourced by the null dust with the energy 
momentum tensor flux and gauge flux mm 

87rGr;;^““^’' = z^[M{v) - 2zq{v)q{v)]6^,s,^, 

= z^q{v)S>^\ ( 10 ) 

where the dot stands for derivative with respect to co¬ 
ordinate V, M{v) and q{v) are the mass and charge of 
a collapsing black brane 


3 Holographic thermalization 

In the previous section, we have got a gravitational col¬ 
lapse solution which describes the collapse of a thin- 
shell of charged dust from the boundary toward the 
bulk interior of asymptotically anti-de Sitter spaces. Ac¬ 
cording to the AdS/CFT correspondence, this process is 
dual to the thermalization of plasma of quarks and glu¬ 
ons which are formed in heavy ion colliders such as the 
RHIC and LHC. To describe the thermalization more 
comprehensive, we should have an initial state and an 
equilibrium state, which is dual to a pure AdS and a 
stationary black brane respectively. Recent investiga¬ 
tions show that as the mass M{v) and charge q{v) are 
written as the smooth functions [Ill[ia[2ll[22] 

-^(^) = ^ + tanh —^ , (11) 

2 V ^ 0 / 

= I f 1 + tanhV (12) 

2 V ^ 0 / 

where vg represents a finite shell thickness, one can con¬ 
struct such a gravity model. In this case, in the limit 
u —> — 00 , the mass vanishes and the backgronnd in 
Eq.® corresponds to a pure AdS space. While in the 
limit u —> 00 , the mass is a constant and the back¬ 
ground corresponds to a phantom Reissner-Nordstrom- 
AdS black brane. Thus for different values of the time 
V, the background in Eq.® stands for different stages 
of gravitational collapse, which represents for different 
stages of thermalization in the dual conformal field the¬ 
ory. 


3.1 Nonlocal probes 

As the gravity model describing the thermalization is 
constructed, we will choose proper observables to probe 


it. Since local observables in the boundary such as ex¬ 
pectation values of the energy-momentum tensor are 
not sensitive to the thermalization process, one needs 
to consider the non-local observables such as the two- 
point correlation functions, expectation values of Wil¬ 
son loops, and entanglement entropy. Recent investiga¬ 
tions show that all of them have similar behavior, thus 
we are interested only in the expectation values of Wil¬ 
son loops. Wilson loop operator is defined as a path or¬ 
dered integral of gauge field over a closed contour, and 
its expectation value is approximated geometrically by 
the AdS/CFT correspondence as [T^[58] 

-^ren (^) 


(W(C)} 


(13) 


where G is the closed contour, U is the minimal bulk 
surface ending on G with Aren its renormalized minimal 
area surface, and a' is the Regge slope parameter. 

Here we are focusing solely on the rectangular space¬ 
like Wilson loop. In this case, the enclosed rectangle can 
always be chosen to be centered at the coordinate origin 
and lying on the Xi — X 2 plane with the assumption 
that the corresponding bulk surface is invariant along 
the X 2 direction. The area of the minimal area surface 
is written usually as |59] 


A = J Vg fiv j iy. 


(14) 


here for notational simplicity, Xi has been replaced by 
X and X2 has been replaced by y. For the Vaidya-like 
AdS black branes in Eq.®, the area of the minimal 
area surface can be written as 

n ^^ ^/l-2z'{x)v'{x)-H{v,z)v'{x)^ 
z{x')'^ 

where I is the boundary septation along x direction. As 
to Eg. dTSl) . we can treat the integral as a Lagrangian C 
and with it we can get the equations of motion of z{x) 
and v{x) 

0 = 4 — 4v'{x)^H{v, z) — 8v'{x)z'{x) — 2z{x)v"{x) 

-\- z{x)v'{x)‘^dzH{v,z), (16) 

0 = v'{x)z'{x)dzH{v, z) + ^v'(x)^dvH(v, z) 

+ v"{x)H{v, z) + z"{x). (17) 

To solve these equations, we need to consider the sym¬ 
metry of the minimal area surface and impose the fol¬ 
lowing boundary conditions 

z(0) = z*, u(0) = u*, u'(0) = z'(0) = 0. (18) 

In addition, the area A is divergent due to its contribu¬ 
tion near the AdS boundary. To eliminate it, we should 
impose a cutoff near the boundary 

= zo,v{^) =to, (19) 
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where zg is the IR radial cut-off and tg is the time that 
the minimal area surface approaches to the boundary, 
which is called as the thermalization time usually. AS 
the divergent part of A is subtracted, the minimal area 
surface, called renormalized minimal area surface, can 
be cast into 

I n 

ri ^2 2 

Aren =2 -. ( 20 ) 

Jg z{xp Zg 

Note that here we have used the conserved equations 
obtained from the Lagrangian C in Ea. dT5|) . and ^ is 
the contribution of the minimal area surface in the AdS 
boundary [IHllIl]. 



Fig. 1 Change of u/T to the charge q. The red line, green 
line and blue line correspond to g = —1, —2, —3 respectively. 


3.2 Numerical results 

In this section, we will solve the equations of motion 
of minimal area surface numerically, and then explore 
how the chemical potential and phantom dark energy 
coefficient affect the thermalization. During the numer¬ 
ics, we will take the shell thickness and UV cut-off as 
vg = 0.01, Zg = 0.01 respectively. 

According to the AdS/CFT correspondence, we know 
that the electromagnetic field in the bulk is dual to the 
chemical potential in the dual quantum field theory, 
therefore we will use the electromagnetic field defined 
in Eq.Q to explore the effect of the chemical poten¬ 
tial on the thermalization process in the AdS boundary. 
When we investigate the effect of the chemical potential 
on the thermalization, the quantity which is physically 
meaningful is the ratio oi u/T, where T is the Hawking 
temperature of the phantom Reissner-Nordstrom-AdS 
black brane. However, as stressed in miiBn], the chem¬ 
ical potential has energy units in the dual field theory 
([m] = 1 /[E]) while defined in Eq.Q is dimension¬ 
less. Thus one has to redefine the electromagnetic field 
as Afe = Afj^/p, where p is a scale with length unit that 
depends on the particular compactification. Therefore, 
through this paper we will use the following ratio 

u ^ ^ 2q{l + r]q‘^) 

T T p(3-?7g2)’ ^ ^ 

to check the effect of the chemical potential on the 
thermalization time, in which we have set the hori¬ 
zon Xh = 1. For simplicity, we will call the ratio in 
Eq. d^Tll as chemical potential hereafter. For the phan¬ 
tom dark energy parameter, we know that it satisfies 
the condition 77 < — 1 , in this paper, we will choose 
rj = —1, —2, —3 in our numerical result. In addition, we 
know that when ry = 1, the solution in Eq.([2]) stands 
for the Reissner-Nordstrom-AdS black brane. To com¬ 
pare conveniently, we will also study the case ry = 1 . 


From Eg. (1211) . we know that for ry = 1, the chemi¬ 
cal potential raises from 0 —>■ 00 provided q changes 
from 0 to v^- While for rj < —1, the chemical poten¬ 
tial is not monotonous as the charge increases, please 
see Figure ([T]). For rj = —1, —2, —3, the chemical poten¬ 
tial increases monotonously for 0 < g < 0.49,0 < g < 
0.35, 0 < g < 0.28 respectively, and monotone decreas¬ 
ing for the other values. To investigate conveniently, 
we will choose g = 0.5,0.7,0.9 in this paper. With 
this choice, from figure © we know that for the neg¬ 
ative phantom dark energy parameters, as the charge 
increases, the chemical potential decreases, which is dif¬ 
ferent from the case rj = 1 where the chemical potential 
increases when the charge increases. 

With the initial conditions in Ea. (fT^ . we will solve 
z{x) and v{x) firstly. When we are interested in the 
effect of chemical potential on the thermalization, the 
phantom dark energy coefficient is fixed, and when we 
are interested in the effect of phantom dark energy 
coefficient on the thermalization, the chemical poten¬ 
tial is fixed. Since different initial time 7 ;* corresponds 
to different stages of the thermalization, we will study 
whether the phantom dark energy coefficient and chem¬ 
ical potential have the same effect on the thermalization 
for different u*. 

The result is listed in Table from it, we can 
observe the following phenomena: 

• For a fixed chemical potential, as the phantom dark 
energy parameter decreases, the thermalization time 
increases, which means that the small dark energy 
parameters delay the thermalization; 

• As the charge increases, the thermalization time in¬ 
creases for the positive dark energy parameter while 
it decreases for the negative parameters. Namely the 
charge has oppositive effect on the thermalization 
time for the positive and negative parameters. How¬ 
ever, as the charge raises, we know that the chem¬ 
ical potential increases for the positive parameter, 
and decreases for the negative parameters. So we 
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= 

-0.888 

r; = 1 
?7 = -1 

77 =-2 

?7 = -3 

q = 0.5 
0.989522 
0.969131 
0.959202 
0.949734 

q = 0.7 
0.999862 
0.9597 
0.94115 
0.923542 

q = 0.9 
1.01398 
0.947474 
0.918376 
0.891606 


77 = 1 

1.34905 

1.37703 

1.41573 

Vi, = 

77 = -1 

1.29452 

1.26975 

1.23855 

-0.555 

77 =-2 

1.26889 

1.22264 

1.16637 


CO 

1 

11 

1.24432 

1.17899 

1.10331 


77 = 1 

1.52816 

1.56579 

1.61914 

Vi, = 

77 = -1 

1.45507 

1.42259 

1.38149 

-0.333 

77 =-2 

1.42132 

1.36125 

1.29001 


CO 

1 

11 

1.38917 

1.30584 

1.21261 


Table 1 The thermalization time to of the minimal area sur¬ 
face probe for different phantom dark energy coefficient and 
different charge at n* = —0.888, —0.555, —0.333 respectively. 

conclude that for all the phantom dark energy pa¬ 
rameter, as the chemical potential raises, the ther¬ 
malization time increases. 

• For different initial thermalization time, the effect of 
the phantom dark energy parameter and the chem¬ 
ical potential on the thermalization are the same. 
Thus we can predict that this effect would be the 
same to the whole thermalization time. 

With the numeric result of z{x), we can further 
get the renormalized minimal area surface defined in 
Eg. d^ . The result is shown in Figure ([H and Figure 
([3|). In Figure ([2|), we fix the chemical potential to inves¬ 
tigate the effect of the phantom dark energy parameter 
on the thermalization, while in Figure Q we fix the 
phantom dark energy parameter to investigate the ef¬ 
fect of chemical potential on the thermalization. We are 
interested in the I independent quantity SA = Aren/l, 
SAq = Aq/l, where Aq is the renormalized minimal area 
surface at the equilibrium state. To observe the ther¬ 
malization conveniently, we plot the quantity 5A — 6Aq. 
In this case, the thermalized state is labeled by the null 
point of 5A — SAq. In each picture, the vertical axis in¬ 
dicates the renormalized minimal area surface while the 
horizontal axis indicates the thermalization time to- 
From Figure (ED, we know that as the phantom 
dark energy parameter decreases, the thermalization 
time increases. That is, the small phantom dark en¬ 
ergy parameter promotes the thermalization. This phe¬ 
nomenon is more obvious as the charge raises. Espe¬ 
cially for the case q = 0.9, we hnd the difference of the 
thermalization curve for different phantom dark energy 
parameters is most obvious. The effect of the chemi¬ 
cal potential on the thermalization is plotted in Figure 
ED- For (a) in Figure ED, background recovers to 
the Vaidya AdS black brane with a chemical poten¬ 
tial which was investigated in [1I1E1E3]- We see that 



(a) q = 0.5 



(b) q = 0.7 



(c) q = 0.9 


Fig. 2 Thermalization of the renormalized minimal area sur¬ 
face for different phantom dark energy parameters with the 
same chemical potential at the same boundary separation 
I = 2.2. The purple line, red line, green line, and blue line 
correspond to tj = 1, —1, —2, —3 respectively. 


as the chemical potential increases, the thermalization 
time enhances. That is, the chemical potential delays 
the thermalization, which is consistent with that ob¬ 
tained in pni221l23] . For (b),(c),(d) in Figure ED, ^e 
find that as the charge increases, the thermalization 
time decreases. This effect is not the same as that in 
(a). Note that for the case rj = — I,—2,—3, the chemical 
potential is monotonously decreasing as the charge in¬ 
creases. Thus we also can conclude that as the chemical 
potential increases, the thermalization time raises too. 
Namely the chemical delays the thermalization always 
though this seems different for the negative parameters 
and positive parameters as the charge increases. 
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(a) ?? = 1 



^0 


(b) 7 ? = -1 



fo 


(c) r? = -2 



(d) 77 = -3 


Fig. 3 Thermalization of the renormalized minimal area sur¬ 
face for different chemical potentials with the same phan¬ 
tom dark energy parameter at the same boundary separation 
I = 2.2. The purple line, red line, and green line correspond 
to 0.5, 0.7, 0.9 respectively. 



Fig. 4 Comparison of the thermalization process between 
the numerical curves and fitting functions. 


d(JA-iA,)/iito 



Fig. 5 Thermalization velocity of the renormalized minimal 
area surface at g = 0.7. The red line, green line, orange line, 
and brown line correspond to 77 = 1, —1, —2, —3 respectively. 


d~(SA-SAp)/dt-„ 



Fig. 6 Thermalization acceleration of the renormalized min¬ 
imal area surface at <7 = 0.7. The red line, green line, orange 
line, and brown line correspond to 77 = 1, —1, —2, —3 respec¬ 
tively. 


In addition, we find the numeric curves for the ther¬ 
malization process can be fitted as analytical functions 
with respect to the thermalization time. Thereafter, we 
mainly focus on the effect of the phantom dark energy 
parameter on the thermalization. As an example, we 
consider the case q = 0.7. The fitting functions of the 
thermalization curves for 77 = 1, —1, —2, —3 can be ex- 
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pressed as 

' -0.532952 + 0.0142505io + 0.27677Ug + 0.111791tg 

- 0.260812t^ + 0.151048ig - 0.0328631tg 
-0.427766 + 0.0156803^0 + 0.00149207ig + 0.442568ig 

- 0.52845t^ + 0.34639215 - 0.093717216 
' -0.380358 + 0.02228361o - 0.183421§ + 0.76406lg 

- 0.9138271^ + 0.6387371^ - 0.1798351§ 
-0.335312 + 0.01267671o - 0.244621lg + 0.73413lg 

- 0.8596831^ + 0.69074215 _ o.219645lg 

Figure (|31) is the comparison of the numerical curves 
and fitting function curves. It is obvious that at the 
order of I®, the thermalization curves can be described 
well by the fitting functions. With the analytical func¬ 
tions of the thermalization curves, we can get the ther¬ 
malization velocity, and thermalization acceleration con¬ 
veniently. From the thermalization velocity curves in 
Figure ([S|), we find the thermalization process is similar 
for different phantom dark energy parameters, namely 
at the middle stage of the thermalization process there 
is a phase transition point which divides the thermal¬ 
ization into an acceleration phase and a deceleration 
phase. This result is reasonable because the thermal¬ 
ization can not accelerate always in order to approach 
to an equilibrium state. The phase transition points can 
also be read off from the null point of the acceleration 
curves, which are plotted in Figure ([HI). For the negative 
phantom dark energy parameter, we find as the param¬ 
eter decreases, the phase transit decreases. For all the 
phantom dark energy parameters, we find the initial 
acceleration of renormalized minimal area surface de¬ 
creases as the state parameter decreases. Correspond¬ 
ingly, the initial thermalization velocity of minimal area 
surface decreases gradually. In addition, from Figure m 
and Figure ([B]) we find the maximum acceleration and 
maximum velocity increase as the phantom dark energy 
parameters decrease. 

4 Conclusions 

We construct a gravitational collapse model which is 
interpolated between a pure AdS space and a phantom 
Reissner-Nordstrom AdS black brane. We probe this 
collapse behavior with the minimal area surface. Ac¬ 
cording to the language of AdS/CFT correspondence, 
this is dual to probe the thermalization behavior of the 
quark gluon plasma by the expectation value of Wil¬ 
son loop in the dual conformal field theory. We focus 
mainly on the effect of phantom dark energy and chem¬ 
ical potential on the thermalization time. We find for a 
fixed chemical potential, the smaller the phantom dark 


energy parameter is, the easier the plasma thermalizes. 
And for a fixed dark energy parameter, the larger the 
chemical potential is, the harder the plasma to thermal- 
ize. For the effect of chemical potential on the thermal¬ 
ization time, there are two points should be stressed: 

(22)» From Figure ([3]), we know that as the charge raises, 
the thermalization time enhances for the case ?7 = 1, 
while it decreases for the case r] = —1,-2,-3. How¬ 
ever, in the dual conformal field theory, the quan¬ 
tity which is meaningful physically is the chemical 
potential not the charge. With the consideration 
that for q > 0.5, the chemical potential raises for 
?7 = 1, while it decreases for r] = —1, —2, —3 as the 
charge grows. Thus in the parameter range which we 
choose, we can conclude that the chemical potential 
delays the thermalization for all the phantom dark 
energy parameters. 

• From Figure CD, we know that the chemical po¬ 
tential is not a monotone function with respect to 
the charge. For 0 < q < 0.49,0 < q < 0.35,0 < 
q < 0.28, the chemical potential is monotonously 
increasing as the charge raises for rj = —1,-2,-3 
respectively. Thus, if we choose the smaller charge, 
e.g. q = 0.2, we find the chemical potential promotes 
the thermalization, which is different from the pre¬ 
vious result that the chemical potential delays the 
thermalization. So strictly speaking, the effect of the 
chemical potential on the thermalization time de¬ 
pends on the dark energy parameters. 

Besides our previous investigation, we also obtain 
the fitting functions of the thermalization curves, and 
with the functions, we further investigate the thermal¬ 
ization velocity and thermalization acceleration. From 
the thermalization velocity curves, we know that at the 
middle stage of the thermalization there is a phase tran¬ 
sition point which separates the thermalization into an 
acceleration phase and a deceleration phase. The phase 
transition point is found to be shifted for different dark 
energy parameters. Especially for r/ = —1,—2,—3, we 
find he phase transit decreases as the phantom dark en¬ 
ergy parameter decreases. Both the maximum acceler¬ 
ation and maximum velocity are found to be increased 
as the dark energy parameters decrease. Though our 
study in this part is not so strict as that in [6T1I62] , 
where the whole thermalization was divided into four 
processes and the analytical function for each process 
was given strictly, our investigation can provide another 
sight to understand the thermalization process. 
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